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We reconstruct the ACDM model for f(T, T) Theory, where T is the torsion scalar and T 
the trace of the energy-momentum tensor. The result shows that the action of ACDM is a 
combination of a linear term, a constant (—2A) and a non-linear term given by the product 

V^Fg [(T 1/3 /16ttG) (167 vGT + T + 8A)j, with F g being a generic function. We show that to main¬ 
tain conservation of energy-momentum tensor should impose that F g [y\ must be linear on the trace 
T. This reconstruction decays in the f(T) Theory for F g = Q, with Q a constant. Our reconstruction 
describes the cosmological eras to the present time. The model present stability within the geomet¬ 
ric and matter perturbations for the choice F g = y, where y = (T 1 ^ 3 /16 ttG) (167tGT + T + 8A), 
except for geometric part to de Sitter model. We impose the first and second laws of thermodynam¬ 
ics to the ACDM and find the condition where they are satisfied, that is, Ta , G e // > 0, however 
where this is not possible for cases where we choose, leading to a breakdown of positive entropy and 
Misner-Sharp energy. 

PACS numbers: 98.80.-k, 95.36.+X, 04.50.Kd 


I. INTRODUCTION 

The description of the gravitational interaction can be done in different ways, the main and best known is that 
using Riemannian geometry as a tool for the formulation of this description. This is known as the Theory of General 
Relativity (GR) of Einstein [If. 

The Riemannian geometry is based on the differential geometry which deals with the so-called space-time as a 
differentiable manifold of dimension four, where we have only the effect of the curvature. There are other geometries 
that generalize this idea. A space-time may have curvature, but also torsion in their structure. There are two 
fundamental concepts in differential geometry. The formulation of this type of geometry has been undertaken and 
gravitation is known as Einstein-Cartan geometry Q. In this theory the gravitational interaction is described by 
both curvature and torsion of space-time, where the torsion is commonly attributed to the inclusion of spin, through 
fractional spin fields. A very particular case of this theory is when we take the identically zero curvature, and then 
only have a space-time with torsion. This type of geometry is known as the Weitzenbock geometry i! , where the 
torsion describes the gravitational interaction. Various analysis can be performed in this type of geometry, which 
is proven to be dynamically equivalent to GR @. In this context, recently a new formulation is started up, and 
generalizes the call Teleparallel Theory (TT), of the space-time where the gravitational interaction is described solely 
by the torsion. 

Through the standard Big Bang Q theory, the ACDM models describe very well the evolution of our universe, in 
a Riemannian geometry within the GR. There are several open issues, but the main one today is the so-called dark 
energy. In order to make the cosmological WMAP data fitting with the theory, it is necessary to introduce an exotic 
component in the equations of GR, i.e, the dark energy. This can be modelled as a perfect fluid with the equation 
of state pde = ljdePde, where the values of lode must be very close to the —1 today. An alternative to that data, 
being consistent with the theory, is the modification of the geometry. A good review can be seen in Q. One of the 
first general possibilities is the well known f(R) Theory [1], where R is the scalar curvature, obtained through the 
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double contraction of the Riemann tensor indices. This theory places the analytic function f(R) in the action, where 
the GR can be reobtained in a certain limit, such as from f(R) = aiR + with <22 —> 0 one gets again GR. This 

theory has proved effective in simulating the evolution of our universe, in various epochs. Other possibilities have 
arisen through the generalisation of the GR, by changing the action. One such change is the case of /(R,T) Theory, 
where T is the trace of the energy-momentum tensor. In this case, the matter content should be taken into account 
as having a kind of interaction with the geometry. 

A direct analogy could be made between theories with only the curvature and the torsion only. As the f(R) theory 
is a generalization of the GR, it is logic to also think of a generalization of the TT, where the analogous object to the 
riemannian scalar curvature, is the torsion scalar T, obtained from contractions between the torsion and contorsion 
tensors. A change in the action of TT is made considering an analytic function f(T) which depends on the torsion 
scalar. That was first thought of a theory arising from the Born-Infcld action [§]. Then, several studies have shown 
the great accordance of this theory with the most varied approaches in Gravitation and Cosmology 0- 

Another recent proposal is to consider not only the torsion scalar in the action, but also the trace of the energy- 
momentum tensor, as an analogy to the /(f?,T) Theory. This theory, called /(T, T) has been formulated recently 
0, still requiring verification of compatibility with the cosmological data and the physical requirements for a good 
cosmological theory. That is the reason to check how should be the functional form of the action of this theory, such 
that the ACDM model is valid. For this we use the scheme the reconstruction method for a modified gravity [121 ]. In 
addition, we make a stability analysis for the ACDM model. 

It also exists a wide interest in studying the thermodynamics of our universe. Various calculus have been done for 
guaranteeing the system to obey the first and second laws of thermodynamics. In the GR, it has been shown that 
the first law of thermodynamics can be shown as dE = TdS + WdV [13|. This law can also be represented in the 
modified versions of gravity, but with an additive content of entropy production, for a non-equilibrium description 
0 . We have done here an analysis which yields the conditions for satisfying the classical thermodynamics laws. We 
adopt the units as kp = c = H = 1 and the Newton constant G~ 1 ^ 2 = Mpi an k = 1.2 x 10 19 GeV fl5[ . 

The paper is organized as follows. In section [TT] we make a brief description of f(T) and /(T, T) theories, with 
the main elements and definitions necessary for their formulation. In section IIIII we have analysed the conservation 
of energy-momentum tensor, which results in strong constraints to the functional form of the action of 
the theory, leading to function /(T, T) with a linear dependence on the trace T. In the IIVI section, we 
use the reconstruction method for the actions to obtain the ACDM model for f(T,T) theory. The result shows 
that the action of ACDM is a combination of a linear term, a constant (— 2A) and a non-linear term 
given by the product y/—TF g [(T 1 / 3 / 167 tG) ( 167 tGT + T + 8A)], with F g being a generic function. In section 
fVl we make the stability analysis of the studied model. The model present stability within the geometric 
and matter perturbations for the choice F g = y, where y = (T 1 / 3 /16 ttG) ( 167 tGT + T + 8A), except for 
geometric part to de Sitter model. We do a thermodynamic analysis for /(T, T) theory in section [Vi] We 
impose the first and second laws of thermodynamics to the ACDM and find the condition where 
they are satisfied, that is, TA,G e ff > 0, however where this is not possible for cases where we choose, 
leading to a breakdown of positive entropy and Misner-Sharp energy. In section I V11 1 we have established 
the conservation of energy-momentum tensor to previous results, showing an inconsistency for this approach. 
We make our final considerations in section IVIIII 


II. f(T) AND f(T, T) GRAVITIES 

In this section we will see the basic preliminary concepts for the reconstruction of theories f(T) and f(T, T) theories 
of gravity. 

In theory /(T), the geometry is determined solely by the matrices that transform the metric of space-time into 
the Minkowski metric. To begin, we define the space-time as a differentiable manifold in which only the torsion is 
non-zero, that is, the curvature is identically zero, then all Riemann tensor components are zero. 

Now we define the line element as 


dS 2 = g gu dx^dx v . (1) 

Taking into account that we can define 1-forms in the co-tangent space of the manifold, and introduce Lorentz 
symmetry in the line element, we can re-write the line element as 

dS 2 = g^dx^dx" = Vab 8 a d b , (2) 

where 9 a = e a g dx^ are 1-forms, index a = 0,...,3 and [ 770 , 5 ] = diag[ 1,-1,—1,-1] is the Minkowski metric. Here the 
Latin indices are related to the co-tangent space and the Greeks ones to the space-time. By the way of writting the 
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line element (J2J , we can establish the following relations g a b = e a ti e b v g^,,, g^ = e a ll e b l/ g a b, e. a ll e a 1 ' = 5”, e a " e b v = 5 b a , 
where e a M is the inverse of matrix called tetrad e a /1 . 

The connection is chosen such that all the Riemann tensor components are identically zero, and one has the 
Weitzenbock connection 0 


V a _ a a 

1 itfi ~ c a • 


( 3 ) 


We can now define a tensor which gives sense of torsion to the space-time 

n, = r % - r v = e a“ - due%) . (4) 

Through the components of the torsion tensor, we can define the contortion tensor components and the tensor S 0 f il ' 


7 y’lMl/ frpflV rpVfl rp {IV\ 

^ a 2 ^ a ^ a ^ a ) i 

ScT = \ + 5» T™ - r a T°* a ) . 

We can also define the analogous object to the scalar curvature in GR, the torsion scalar 

rp rpOL Q {IV ^_rp{lV(T rp , ^ rp{Ll/<Jrp _ rp (JflrpV 

1 — 1 {iv D ot — ^ 1 {iva T ^ 1 av{i 1 G 1 v{i ’ 


( 5 ) 

( 6 ) 

( 7 ) 


It is this object that plays the curvature scalar role in GR, and that should also form the action of f(T) theory. The 
action of f(T) theory is constructed in a way that we have a linear term in the torsion scalar, another containing the 
correction term to the TT and another term related to the material content. Then we write the action as 


S= 1^G / ^ xe [ T + /( T ) + £ ™]’- 


( 8 ) 


with e = det[e a ] = = y/detlg^], G the Newton’s constant, and setting c to unity. In f(T) theory, tetrads 

are dynamic fields, then doing the functional variation of the action 0 in relation to them, one gets the following 
equations of motion 


(1 + / T ) \e-%(ee a °S a ^) e a x T^ x S a ^] + e a °S^d.Tf TT + -e/[T + /] = 4 ttG<© 0 


( 9 ) 


where we use the nomenclature fr = df/dT and /tt = d 2 f/dT 2 , Q a v represents the components of the matter 
energy-momentum tensor. 

Let’s take the example of the Friedmann-Lemaitre-Robertson-Walker (FLRW) universe with flat spatial section 

ds 2 = dt 2 — a 2 (t) ( dx 2 + dy 2 + dz 2 ) , ( 10 ) 

where a(t) is the scale factor. The Hubble parameter is given by H(t) = a(t)/a(t). Now we specify our choice of 
tetrads as 


K/J = diag[l,a{t), a(t), a(t)]. ( 11 ) 

With this choice, we can represent the line element (121) through a set of 1-forms [0° = dt , 6 1 = a(t)dx,d 2 = 
a(t)dy,6 3 = a{t)dz\. Thus, the equations of motion © for f(T) theory, taking the material content as a perfect 
fluid Q" = diag[p mt , - p m t, ~Pmt, ~Pmt ], are given as follows 

H 2 = — 77 — Prat 2H 2 f T (12) 

3 6 

• 4:7TG(p mt T Pmt ) 

-- 1 + f T -l2H 2 f T T [ j 

The subscribed p m t and p m t means the density and pressure of the total matter in the universe, we consider here only 
the components of the baryonic matter as {p m ,p m } and radiation as {p r ,p r }. 

Here, the torsion scalar is obtained by the definitions 0-0 with m, resulting in 

T = -6 H 2 , 


(14) 
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We now see clearly that the equations of motion of f(T) theory are identical to that of the GR, and the equations of 
Friedmann (flat spatial section), when the nonlinear term are zero, i.e making f(T) = f T = frr = 0 in (fl2l) - (fl3l) . 

Now we can present the most recent generalization of the f(T) theory. Following the analogy of the generalization 
of the f(R) theory to /(i?,T), where T is the trace of the energy-momentum tensor. Here, we can also introduce 
in the action an analytic function that depends not only on the torsion scalar, but also on the trace T. For such a 
theory, the action is then given by 

5 = TF~r< l d * x e t T + T ) + l ^GC m ] , (15) 

lb 7T Gr J 

where f(T,T) is an arbitrary analytical function of the torsion scalar T and of the trace T of the matter energy- 
momentum tensor 0^, and C m is the matter Lagrangian density. 

Here we must consider that the Lagrangian density C m depends only on tetrads and not on its derivatives. For the 
energy-momentum tensor of a perfect fluid we have the following trace 


T = 0/ = Pmt ~ 3pmt • (16) 

We can then make the functional variation of the action m with respect to the tetrads, resulting in the following 
equations of motion El 


(1 + f T ) [e~%(ee a a S a ^) - e^T^S™] + (frrd.T + fr+d^T) e a a S a ^ 


+e„ 


f + T 


-h 


*© 
v -' n 


‘ Pm& a 


= AnGe “0„ 


(17) 


where fj- = df /dT and frr = d 2 f/dTdT . Here it is evident that the particular case where the function / depend 
only on the torsion scalar T, i.e, / = /(T)[3S|, the equation (fTTl) reduces to the equation of motion of f(T) theory in 

©• 

Taking again the matter content as a perfect fluid and choosing the diagonal tetrads in m, the equations of motion 
of f(T, T) theory for a flat FLRW universe, are given by 

3 H 2 = 8 nGpmt - 1 (/ + 12 H 2 f T ) + fr (Pmt + Pmt) , (18) 

H = -47 tG (p mt T Pmt) R (fr - 12 H 2 f TT ) - H ( p mt - 3 p mt ) f TT - fr ^ Pmt + Pmt ^j . (19) 

In this paper, in the next section, we propose to reconstruct the ACDM model and study the stability of the de Sitter 
and power law solutions. 


III. CONSERVATION LAWS TO /(T, T) THEORY 


This section is devoted to establishing a conservation law for the f(T,T) theory. For this, we must describe the 
equations of motion in a covariant form. Let us multiply the equation of motion 03 by PnaG’ij (to have the indexes p 
and v the free in the end) and using the identity g^e 0 ll [e~ 1 d\{ee a a S a aX ) - e a a T x rya S x ' ia } = (1/2) [G^ - (1/2 )g flv T\, 
where G^ v is the Einstein tensor, we can rewrite our equation of motion as 

i(l + f T ) (g^ - \g^ + S v x (,frrdxT + frrdxT) + \g^ v (f + T)- X -f T (0„„ + g^p mt ) = 4t rG0^ . (20) 
Now we take the divergence all two sides and isolate taking into account that V= 0, what gives us 

V M 0/ = (47rG+(l/2)/ r ) { {fTTdaT + f Trd ^ e 0 „[e" 1 aA(ee 0 “5 a ^) - e a a T\ a S^] -1(1 + f T )d v T 
+V m 5/ a (frrdxT + frrdxT) + S^ x (j frrrd^Td x T + frrrd^TdxT + frrT^dxT + frrrd^TdxT 
+frrrd„TdxT + /ttV^aT) + J (f r d v T + f r d v T + d v T) - l (frrdxT + frrdxT) (0/ + 8(f Pmt ) 


.j f'T^vVmt 


( 21 ) 
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This equation gives us two constraints, one for v = 0 and other for v = 1, 2,3, as follows 
Pmt A {^Pvnt T Pmt ) = 


dT dT 

2(p m t T Pmt') ( frT-jj- + 


, „ dp m t \ 

-/r^ + 2/T- 


dt j’ 


0 = 


Pmt r, dT , r dT 

— fTT^ + frr^ 


The second constraint results in 


dt 


[It] = 0, 


( 22 ) 

(23) 


(24) 


fixing /(T, T) as a linear function of the trace T or a constant, which is in agreement with the results obtained recently 
in [32[ for the analogous f(R, T) Theory. Already through the first constraint, if we are to maintain the conservation 
of energy-momentum tensor, we must impose p m t + 3 H(p mt + p m t) = 0, as usually obtained in the ACDM model. 
This provides us 


2 (P mt “t" Pmt) 


f dT +f dT 

/TT ^ + /rr df 


- /r 


dr 

dt 


2 fr 


dpmt 

dt 


= 0. 


(25) 


When we do the same analysis in f(R, T) Gravity, we also have a restriction on the functional form of the action 

0. 


We will guard these results and reconstruct the ACDM model taking into account these two constraint. 


IV. RECONSTRCTION OF ACDM MODEL IN f(T,T) THEORY 


In this section we use the reconstruction method, through a particular model, for obtaining what should be the 
functional form of the function f(T, T). This method basically consist in choosing a model consistent with the 
cosmological data, and use the imposition of the characteristic equation for this specific model, must be satisfied at 
any time of the evolution of our universe. So with an imposition, we integrate the equation of motion in order to make 
the model being valid. This results in a functional form of fixed function that makes up the action of the theory. The 
use of this method in some cases of modified gravity can be seen in BUS- 

NOW, for reconstructing the ACDM model, it is necessary to impose an equation of motion of this model 

3 H 2 = 8nGp mt + A , (26) 


where A is the cosmological constant. In this model, the matter is described through a perfect fluid formulation, 
where the pressure satisfies to the following equation of state p m t = to m p m + u r p r = (1/3 )p r , with p m and p r being 
the matter and the radiation densities. The trace of the energy-momentum tensor m is 


r — Pmt 3Pmt — (1 dca m )p m T (1 3uj r )p r — Pm • 


(27) 


By inverting the relation (1271) for the energy matter density in terms of the trace and using (|14D . imposing the 
ACDM model from the equation (HI Hi we can substitute in the equation of motion (fl8l) of the /(T, T) theory, 
yielding 


A = “(/ — 2T/ t ) 



T + 


1 

2?rG 




Integrating this differential equation we obtain the following action function 


f(T, T) = —2A + V—TFg 


- J,l/ 3 

— (16.GT+T + 8A) , 


(28) 


(29) 


where F g [x] is a generic arbitrary function of its argument x. Taking into account the constraint (1241) . we see that 
T'gfa;] = x , linear function in x. The function F g [x\ can also take a constant value. We see that the reconstruction of 
the action (fl5l) is given by 


S = 


1 


167rG 
167tG . 


J d 4 xe [T + f(T, T) + 167rG£ m ] 

yl/3 


d 4 a 


T-2A + V—T 


107rG 


(16ttGT + T + 8A) 


+ 16ttG£„ 


(30) 
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We can now see that this action generalizes the f(T) theory, because being able to have terms representing the 
interaction between the torsion and matter, generated by the function F g in (1^1) . We also see clearly that this 
generalization is not for any functional form between T and 7”, but only by the product given in the third term of 
the above action. This restricts a lot the functional form of the theory, when it is imposes the validity of the ACDM 
model, as we have here. Another important observation is that this action lies_in the f(T) when we properly choose 
F g [(T 1 / 3 /167rG) (167 tGT + T + 8A)] = Q, where Q is a constant given in [l8j |. In the next section, we will use a 
model for this generic function, where F g [y\ = y , with y = (T 1/,3 /167tG) (167tGT + T + 8A). 

Here it is still some important observations. Our model is able to reproduce well the eras that the evolution of 
our universe goes. For exemple, with the imposition (1261) we see daily that the era of radiation, where we write now 
3 H 2 w 87r Gp r + A, and of the matter, where 3 H 2 ft 87 t Gp m + A, are now well reproduced. In radiation era, as 
the term of the radiation density dominates over the matter term, the generic function is being approximated by 
F g [(T 1 / 3 /167tG) (T + 8A)]. Already in the matter era, as now the dominant term is the matter density, we have the 
same functional form of the action (1301) . In the dark energy era, it’s approach 3 H 2 ft A, or equivalently —(T/2) ft! A, 
Which brings us to the same functional form of radiation era for F g , but now T ft —2A. Here inflation can be 
described similarly, whereas the component that dominates is the vacuum energy density. 

As seen, our reconstruction is compatible with the most varied eras of cosmological evolution of 
our universe, leading us to believe that this model should also be compatible with the cosmological 
experimental data. We know at least that the acceleration of the universe at this stage is consistent 
with our model, for the very reconstruction is done to satisfy this condition. And we still have 
compatibility with the measures of type IA SN, where ujde ~ —1 |36|| . We can at least see that there 
is a compatibility with ultra-stiff fluid where Lo m = 1, made in a recent analysis |35j. 

One last important observation is that for the reconstruction satisfy the first constraint (l?5l) . we have to fix the 
trace in terms of torsion scalar as follows 

T = —f(11 + 4608GV) T — 8Al . (31) 

167tG 

This shows us that the validity of the reconstruction of the ACDM model is subject to a fixation of the energy- 
momentum tensor trace, given in (1311) . This show us that the /(T, T) theory should always be considered a theory 
in which the trace of energy-momentum tensor depends linearly on the torsion scalar. However, this does not permit 
us to conclude that the f(T, T) Gravity must always fall back on /(T) Gravity, because for what /(T,T) fall back in 
/(T), it is necessary that F g = Q in (EH1) . and does not have a linear dependence on the trace, as we have here. 

It may be possible to escape from the result that the action should have linear dependence on the 
trace T, which is a direct implication of the conservation of the energy-momentum tensor given in 
(124)) . We can think about the possibility of non-conservation of energy-momentum tensor, arising from 
the non-minimal coupling between matter and torsion, similar to the case of theories modified only 
curvature j33-1, 1341 ]. This should be interpreted as a possibility that the movement of the material content 
is not more geodesic having the appearance of an extra force term. This can also come to modify the 
bending of light and add a term in acceleration exerted by gravity, astrophysical implications may 
have to help in models of dark matter. This approach goes beyond our analysis here and should be 
done in future work. 

In the next section, we will show the stability of the ACDM model reconstructed here. 


V. STABILITY OF DE SITTER AND POWER-LAW SOLUTIONS 

In order to establish the validity of the ACDM model, let us make a simple stability analysis. To do so, we will do 
a small perturbation in the geometry and the matter 

H(t) = h(t)[ 1 + 5(f)] , pmt{t ) = pmth(t)[ 1 + S m (t)] , (32) 

where {h(t), p m th{t)} is an exact solution of the equations of motion for call background (obey to (fI71) - (fT51) l and 
S(t), 5 m {t) « 1. With this, the torsion scalar is given by 

T = -6H 2 = —6h 2 (l + S) 2 = T 0 (f)[l + 2 S(t)} , (33) 

where we have To(t) = — 6h 2 (t). The trace of the energy-momentum tensor is given by 


T(f) = Pmt{t) - 3 Pmt{t) = (1 - 3uj m )p m {t) + (1 - 3u) r )p r (t) = p m (t ) = T 0 {t)[l + 6 m (t)] , 


(34) 
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where 7o (t) = p m h(t)- The perturbations of the torsion scalar and the trace of the energy-momentum tensor are 

5T(t) =T-T 0 = 2T 0 (t)8(t) , 5T{t) = T- To = To(t)5 m (t) . (35) 

Let us express the function /(T, T) in the Taylor series until the first order around the point [T 0 (f), 7o(t)], by 


/(T,T)~/(T 0 ,To) + 


its derivatives, that are obtained deriving (j36l) . can be expressed as 


[df] 

ST + 

'df 

[d'T \ 

t=t 0 

[dT \ 


ST. 


Jr=r 0 


fATT) « /t(T 0 ,To) + 
fr(T,T) » MT 0 ,To) + 


d 2 f 


dT 2 


ST - 


JT=T 0 

d 2 f 


d 2 f 

dTdT 

St 


ST , 


Jt=t 0 ,t=t 0 

d 2 f 


Jt=To,T=To 


<9T 2 


ST. 


JT=To 


dTdT_ 

By substituting these expressions into ED- collecting the term until the first order, we get the following result 

3/i 2 (l + 28) = 8nGp mth (l + S m ) - [/ 0 + f To ST + f To St] + I2h 2 (l + 28) x 

X [ho + hoT 0 ST + hoT 0 ST] + | + [ho + hoToST + fr 0 T 0 ST] (^Pmh + ^Prh^j (1 + Sm) 
that replacing p m h = To, = —T 0 /6 and p r h = — [(1/87tG)((To/2) + A) + To] becomes 


(36) 

(37) 

(38) 


(39) 


-To < 1 + fr 0 + 2T 0 /t 0 t 0 - 


To 


1 (T 


2t tG \ 2 


-A 


hoTo D S(t) = < Y + A 


5To + ^(t +a 


h 0 + To%hoTo - gT 


T 0 + 


l T 


27 tG V 2 


+ A 


h 0 To U m (t). 


(40) 


We maintain for the moment this result. Now we have to express the perturbation S(t) in terms of 8 m {t ). To do so, 
let us take the perturbation of the conservation equation of the energy-momentum tensor 


Pmt T 3 H ( Pmt T Vmt) — 0 

Pmth (1 + 8 m ) + pmthSm + 3/l(l + 5) ( p mh H - Prh ) (1 ”1“ ^m) 0 


m = -. 


Pmth 


-Sm(t) = -- 




[(To/2) + A] 


-8 m {t) ■ 


3 h (Pmth + Pmth) 8nG\/—T q { (l/8irG) [(To/2) + A] + 7o} 

Here, the “dot” denotes the derivative with respect to time t. Now, by inserting © in ED- one gets 

Fi(t) 


Sm = exp 


f '<*> = I + a- 6 


F 2 (t) 
1 


dt 


5T 0 


1 T 


7tG V 2 


A 


ho + ToTihoTo - g 7 o 


To 


1 (Tr 


2irG \ 2 


+ A 


F 2 (t) = \/— T 0 (^V6To + 2A^j {— 67tG [1 -I- ho + 2T 0 /t o t o ] + (47tG7o + T 0 + 2A) h 0 T 0 } j 
x (4t tGT 0 + T 0 + 2A) 


floTo - 

24?rG x 


(41) 

(42) 

(43) 

(44) 

(45) 


This equation is valid for specific models, the ones we have chosen, the de Sitter and power-law models. We will 
substitute the characteristic of these models in the above equation. 

Let us first analyse the stability of de Sitter solution, where the characteristic is given by h = ho, yielding through 
the equation of continuity 


Pmhit) = po exp[—3/i 0 f] . 


(46) 









































The torsion scalar is constant and is given by Tq = —6/ijj. The trace of the energy-momentum tensor is given by 

% = po exp[—3/i 0 i]. (47) 

Now we have to use the reconstruction of f(T,T) for ACDM in (1^1) . specifying the case of de Sitter solution. We 
first consider the time depending torsion scalar, i.e, variable torsion scalar, and then fix for the de Sitter solution, we 
can then correctly calculate the derivatives of the function /(T, T). 

For the de Sitter solution, we take one case for action (l30l) . We make the case where F g [y\ = y, with y = 

(T 0 1/3 /16ttG) (167tGT 0 + To + 8A). 

For this case, we have the following integral in C3l) 


S m (t) = exp i - 


4 x 6 1 / 3 (/ig /3 + 6 1 / 6 7rG)(3/ig 
72/37777177 7" 


- A) 


/q/,2 


/ig /3 (6 5/6 /jQ — 67rG/ig' /3 )(3v / 6/io ~ A) 


and from m 


(48) 


4 x 6 1 / 3 (fep /3 + 6 1 / 6 7rG)(3feg - A) [3 Kj - A - 27rGpoe~ 3,tot ] 
h b J 3 (6 5 / e h,Q — 6nGh^ 3 )(3\/6hQ — A) 


(49) 


We do a numerical graph representing the temporal evolution of these perturbations in Figure [T| We can see 
that the perturbations rapidly decrease to zero, showing a possible stability for the solution of de Sitter. There is a 
serious problem here. The graphical representation to S(t) in figure [I] shows that despite this function quickly drop 
to zero, the initial values are of the order to 10 41 , thus contradicting the initial assumption for perturbation which is 
S(t) « 1. This shows us an impossibility for stability to the model of de Sitter, due precisely to the initial values for 
perturbation of the geometric part. 




Figure ll The left graph shows the temporal evolution of 5 m (t) (in red) and the of the right the temporal evolution of <5 m (f) (in blue) both for 
the case F g [y] — y. with y = (t) / :1 /1 OttC') (I GttG'To T To + 8A), the two for the particular case of de Sitter solution. The parameters are chosen 
as {N = 20,h 0 = 2.1 x 0.7 X lO -42 ,^ 5 = 1, p 0 = 0.1 X 10 _121 ,G = (1.2) 2 x 10 _38 ,A = 10 -42 }. 


Doing the same procedure for the case of power law, where now h, 0 (t) = a/1, the equation (l43l) . for the case F g =y 
with y = (Tg 1 ^ 3 /167 tG) (167tG7o + T 0 + 8A), becomes 


a) 


pl2 

3a 


5 m {t) = exp | J y^dt 
Fpi\ = -4 x 6 1/3 

F pl2 =a 2 (3V6a 2 -At 2 ) 
From (l42l) . with (l50l) . we find 


2 + 6 ll6 irGt 2 (|) 1/3 
6 5 / 6 a 2 -6nGt 2 (|) V3 


(At 2 — 3a 2 ) 


(At 2 - 3a 2 ) 


3a; 


ZirGpot 2 


(50) 

(51) 

(52) 


m 


Fpl2 


at 


At 2 — 3 a 2 

¥+^G P0 (±y 3a ^j 


(53) 
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We numerically represent the evolution of perturbations 6 m (t) and S(t) in the figure [2] Again we can see that the 
perturbations quickly decay to zero, showing the stability of this solution. 




Figure 2: The left graph shows the temporal evolution of S m (t) (in red) and the of the right for temporal evolution of S(t) (in blue) both for 
the case E g [y] — y, with y = (T ( j / :i /1 6ttG) ( IGrrGTo + Tq + 8A), the two for the particular case of power law solution. The parameters are chosen 
as {a = 2, to = 3, h 0 = 2.1 x 0.7 x Id -42 , <5™ = l,p 0 = 0.1 X 10 -121 , G = (1.2) 2 x 10 -38 , A = 10” 42 }. 


In the next section, we will study the two laws of thermodynamics for the ACDM model. 


VI. THERMODYNAMIC LAWS ACCORDING TO ACDM MODEL 

By analogy to the thermodynamics of black holes, the thermodynamics of the cosmological models can be realised, 
verifying if these models satisfy the classical thermodynamics. This has been done in the context of GR |2cj. It has 
also been of wide interest the study of these laws for the case of the theories of modified gravities [21] . In the context 
of modified gravity it is common to use a description of non-equilibrium thermodynamics for representing the physical 
system [22j . we will use this representation here. 

In this section we will investigate the conditions of satisfying the first and second laws of thermodynamics. 


A. First law 

We will redefine now the action of the theory /(T, T) as 

S = j d A xe [. f(T , T) + 167rGX m ]. (54) 

This is necessary to get the possibility for defining an effective Newton constant G e //. Now the equations of motion 
are rewritten as 


3H 2 = 8nG eff (p mt +p DE ) , (55) 

H = -AirG e ff (p m t + Pmt + Pde + Pde) 1 (56) 


with 


Geff = 


Pde = 


G_ 

It 


1 + 


/r 


167tG 


1 


16 -rrGefffT 
Pde = Pmt + Pmt ~ Pde — 


fpPmt rj / 


4?r Geffh L 


12 H 2 Hf TT - IttH ( Pmt - 3 Pmt) 


(57) 

(58) 

(59) 


We impose the conservation of the matter sector p mt + 3 H(p mt +p m t ) = 0, as discussed in the sections Hill and HV1 
But as we are doing a non-equilibrium thermodynamics description, the dark energy sector is does not conserve, 


Pde + 3 H(pde +Pde ) 


3 H 2 d / 1 \ 
8tt dt \Geff) ’ 


(60) 
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where we used (I55l) - (l59|) . This concords with the theory analogue to f(T, T), in [23|, and also for the particular case 
where / = f(T) in m We still should compare our results with the ones of f(R) theory, in [24]). It is also clear here 
that, when G e // = G (in the particular case f(T, T) = T) the dark energy sector is conserved, the TT is recovered, 
or the GR analogously, within an equilibrium description of thermodynamics. 

Now, we establish the basic tools for the first law of thermodynamics. The line element m can be rewritten as 

dS 2 = h^dx^dx" + r 2 (dd 2 + sin 2 ddq i 2 ) , (61) 


where we have the metric of 2-dimensional space [h^ = [1,— a 2 (t)], for p, v = 0,1 and x 1 = t, x 2 = r, a new 
radial coordinate r = ra(f), with r being the usual radial coordinate r 2 = x 2 + y 2 + z 2 . Here the space-time can be 
decomposed in a 2-dimensional space with the metric and other 2-dimensional space with is a 2-sphere, with the 
line element dtt 2 = dd 2 + sin 2 ddcj) 2 . 


Through this line element (15TT) . we can calculate the apparent horizon 41] and the associated temperature to this. 
The apparent horizon is obtained from the expression h ,,v dS A d v r A = 0, which results in 


1 

Ta = h- 


(62) 


The temperature is calculated by the expression 


Ta = 


1 


47 x\J—h 


yf—hh^duTA 


(63) 


where h = det[h llv \. The temperature is given by 


T a = 


r_A_ 

47r 


(H + 2 H 2 ) , 


(64) 


where, from now, we impose H + 2 H 2 > 0 for getting a definite positive temperature. 

Using HMD and m, we can calculate the derivative of ?a 

^ = 4 tt G e ffr\ (ptot + Ptot ) , (65) 

with ptot = Pmt + Pde and p to t = Pmt+PDE- Now, we will take the entropy related to apparent horizon, for modified 
gravity theories [25], that is Sa = A/{A G e ff), with A = Awr 2 ^. Here, appears a new restriction for the thermodynamic 
system, G e // > 0 for getting a definite positive entropy Sa > 0. Considering this entropy, and using G e // in (l57l) 
and (1551) . we get the following differential 


dSA = 2 -kta 


Atvta (Ptot + ptot) dt + -r A d 


G , 


eff n 


Now we take the Misner-Sharp energy for modified gravities Ems = rA/(2G e //) = V p t ot\M 
where, using (|65|) . the differential gives rise to 


dE MS = 2nr 2 A (p tot + p tot ) dt + -r A d 


G, 


eff 


Using ([62]) in (El, considering El, we can rewrite the temperature as 


T a = 


t_a_ 

47T 


1 

dt \ta 


_ 1 / _ 1 dr A \ 


2nrA 


2 dt J 


Making the product (1551) for (1551) we have 


( 66 ) 

|, with V = (4/3)tt r 3 A , 

(67) 

( 68 ) 


TaASa = —dE M s + 27t r\ ( ptot - Ptot ) dr a + 67rr^ ( p to t + Ptot ) dt - ATT? 2 A ptd?A 
+ -jVA (1 + 2e?aTa) d 


G, 


eff 


(69) 
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Defining the work density given as W = (1/2 ){p to t — Ptat) Hi using ([55]) and (l65l) . we can rewrite (l69l) as 

T A dS A + T A dS p = -dE MS + WdV , (70) 

T A dS p = ~^r A (1 + 2tt? a T a ) d (^~) - ( 71 ) 

where the term T A dS p can be interpreted as coming from an entropy production. Considering T A > 0, we eliminate 
the singularity in the entropy production in H = —2 H 2 (T A = 0). Here, we see that the result concords with the 
particular case where / = /(T), in 0- Also in the analogue theory f(R,T) in jUJ. This entropy production dS p 
is directly linked to the non-conservation of the dark sector in (ICJUl) . When we fix the particular case f(T, T) = T 
( G e ff = G), we recover a theory without entropy production, i.e, the TT. 

We take the effective constant (l57l) to the power law solution, which provides the 


Geff(t) 


/ \ 3a 

eG (i) 

8 x 6 1,/6 7rGt 2 (f) 1/3 — 3a 2 


2 « 2 ( A (i) 

3a \ / \ 3a 

+ 2ttGpqJ — 33a 2 


(72) 


We represent the temporal evolution of the effective Newton constant G e //(f) in the figure [5] we show what for the 
case of power law, G e // assumes negative values, precluding a definition of Misner-Sharp energy and entropy positives. 



0 1 x 10 21 2 x 10 21 3 x 10 21 4 x 10 21 5xl0 21 6xl0 21 


t [ 1/GeV] 


Figure 31 The graph shows the temporal evolution of G e ff(t) for the case F g [y] — y, with y — (t/^ 3 /167tG) (167tG7o + To + 8A), for the 

particular case of power law solution. The parameters are chosen as 
{a = 2, to = 3, h 0 = 2.1 x 0.7 X lO -42 ,^ = l,p a = 0.1 X 10” 121 , G = (1.2) 2 x 10 -38 , A = 10 -42 }. 

So we can conclude that this type of behavior suggests (it we have not tested all possibilities of cases) that the 
/(T, T) Gravity, does not have a well-established thermodynamic with all the thermodynamic quantities physical 
well-defined every time interval. 


B. Second law 

We will impose the second law of thermodynamics to f(T,T). To do so, let us consider the Gibbs equation in terms 
of matter fluid and total energy given by 


TtotdStot (Ptot T Ptot ) dV T Vdpt 0 t , 


(73) 


where T to t and Stot are the temperature and total entropy associated to the system, respectively. According to the 
second law of thermodynamics, we must have 


dStot ^ dS A ^ dS p ^ q 
dt dt dt ~ 

Knowing that Ems = V Ptot and using CoD and (1731) we can write the left side of (1771) as 


(74) 


dStot 


dS A 

dt 


dS p 

dt 


1 

T A Ttot 


(Pfof 4“ Ptot ) 


dV 

dt 



Ttot 

— 


+ V(T a - T tot ) 
at 


dt 


(75) 
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For Ta = T to t, using (156)) and dV = Airr\drA we get, for (1751) . 

dS tot dS A dS p _ 1 / dr A \ 2 

dt + dt + dt 2T A G e ff V dt ) ’ [ } 

with G e ff > 0, as previously done. The observation is that the temperature of the apparent horizon in (1641) . depends 
on the components of the universe, through (155]) and (1561) . including the photons of CMB with 2.73K. Assuming 
that the total temperature of the system Ttot is equal to the one of the apparent horizon, i.e, this temperature also 
includes the photons of CMB. 

Here are the expressions of the second law similar to particular cases, or analogs f(T) and f(R, T), but we are 
still restricted to the problem of the previous section, in which we can not set G e // always positive, which leads to a 
possibility of negative entropy, which is physically inconsistent for the cases studied here. 


VII. COMPARISON WITH PREVIOUS RESULTS 


We noticed a great need to compare our results with the first work in this extension of the f(T) Gravity. The 
observation that we have done here is that the first construction of f(T,T) theory has been done by Kiani and Nozari 
[29l ]. But do to the ambiguity of the choice of the Lagrangian density for the perfect fluid, with the Lagrangian density 
being C m = p (our choice) or C m = —p (see ®), they use the last description. This conduces to a new approach 
possibly equivalent to the same theory. They make an analysis of the stability for the de Sitter model, but as they 
did not reconstruct the action, they considered the model as a generic one /(T, T) = k\T + k- 2 T m T n . 

Here we will verify the consistency of this approach to the conservation of energy-momentum tensor. Tak¬ 
ing the equation of motion (11) in [29], considering (13), we multiply by e“„ and considering the identity 

= ( 1 / 2 ) [G„" - ( 1 / 2 )^T], we have 

S^frrdoT + \f T (g/ - 1%t\ + K [\f + \pmtfr ) = (^G + ^/ r ) 0/ . (77) 

Now we take the total divergence of the above equation of motion, and isolating V,,0/, taking into account that 
V^Gj/ = 0 , we have the following result 


(47rG + (5/4)/ r ) ^ VfJ- 


[/ttV^T + d a T (frrrdpT + 


+ 0 (/tt^T + /tt^T) ( G/ - -6ZT) + - f r d v T + - [frd.pmt + Pmt (. hrd v T + frrd.T) 
(f TT d^T + .frrd^T) 0/ 1. 


(78) 


Then arise two constraint here, one for v = 0 and outher for ^ = 1,2, 3, 

Pmt + 3H(p mt + Pmt ) = 4 I [pmtfTT + 12(1 — 2 H 2 )/ T t\ + [PmtfrT + /t + 12(1 — H 2 )f Tr ] — + 6frPmt}j?ty) 

0 = 3(77 2 - a) Stt—tt + /tt-tt • ( 80 ) 


The second constraint is satisfied, in general, to ( d/dt)fr = 0, fixing function f(T,T) as a linear function on the 
torsion scalar T. This returns to a trace coupling with Teleparallel Theory, which is not the intention here. Thus, 
the description of the theory as shown in [29(, is not consistent with the conservation of energy-momentum tensor 
presented here. We must then discard here that description by the density of the fluid. 

An important observation is that this result does not mean that the description of Lagrangian fluid through the 
density is incorrect, but that the specific choice of the set of diagonal tetrad, [e b ] = [ 1 , a(t), a(t), a(i)], is that it can 
not provide a consistent result. It should be investigated in another future work, as well as consistency between two 
description possibilities, one by density and the other by fluid pressure. 
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VIII. CONCLUSION 

The f(T, T) theory can be a good opportunity to test whether we can describe the dark energy by changing the 
TT equations of motion. To do this, it is first necessary to test the most common models of modern cosmology, the 
famous ACDM one for exemple. 

We reconstruct the gravitational action of this model and have done a brief analysis of their stability. The result is 
resumed as follows. 

The reconstruction of ACDM model gives rise to an action as combination of a linear term in the torsion scalar, 
a constant (—2A) and another, which is the interaction between matter and torsion, described by the product 
\J—TF g [(T 1/,3 /167 tG) (167tGT + T + 8A)1. This reconstruction decays in the f(T) theory, by appropriately choosing 
the F g = Q, with Q a constant given in [18}. This reconstruction severely restricted the functional form of action, 
but well reproduces the cosmological eras. 

We do the analysis of the energy-momentum tensor conservation through the obtained reconstruction 
model, resulting in the function F g [x] must be linear in its argument, so that conservation is satisfied as 
the usual ACDM model. With this action should not contain terms of higher orders and is restricted 
to linearity in the trace T. This model does not prevent accordance with observational data, and for 
the particular case of ultra-stiff fluid, where u> m = 1 (see [35}). 

We have done the analysis of the stability of this model, for a particular choice F g = y , and we found that the 
model is stable for de Sitter and power law solutions, except for geometric part to de Sitter one. 

We finish analysing the thermodynamics for the ACDM model. We find that the imposition of the thermodynamic 
laws are satisfied for the conditions TA, G e // > 0, which can not be met for the analysis here. We found an impossibility 
to define positive entropy and Misner-Sharp energy in any interval of physical time, as well as prevent ghost structure 
in the quantized theory [28}. This suggests that the cosmological thermodynamic theory to /(T, T) may not be well 
established fundamentally. 

Our perspectives are listed as follows. 

We use here the reconstruction method to cosmological models, which can then be examined re¬ 
construction to other known models like Holographic and Chaplygin Gas for example. These new 
approaches can not modify the main result of the energy-momentum tensor conservation, the action is 
linear in the trace T, more should cooperate with the acceleration of the universe at this stage. This 
will be covered in detail in a next work. 

Another possibility approach is when we consider the space-time connection, regardless of the set 
of tetrads, then one of the dynamic fields. That would be a new version of /(T, T) Cosmology in the 
Palatini’s formalism, which may reveal new results in the reconstruction of the ACDM model. This 
approach can also reveal a different constraint to the conservation of energy tensor-momentum. Using 
the formalism reconstruction from the LQC model, as seen in [37}, we can generalize to a non-minimum 
coupling with the matter. We see in [371 ] the reconstruction by an inverse method considering also 
p = p(lZ) or yet 1Z = TZ{T), what is obtained in terms of the torsion scalar on (l3lTl . Could still think 
of the non-conservation of momentum-energy tensor as a consequence of quantum gravity anomaly 
[? ], then giving the possibility of action contain non-linear terms on the trace T. This should be 
addressed in future work. The non conservation may also be arising from quantum effects such as 
creating particles in the expansion of the universe, which also occurs in Rastall’s Cosmology [38|] . 
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